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SUMMARY 


An unsteady potential flow analysis, which accounts for the effects of blade 
geometry and steady turning, is being developed to predict aerodynamic forces and 
moments associated with free-vibration or flutter phenomena in the fan, compressor, 
or turbine stages of modern jet engines. Based on the assumption of small-amplitude 
blade motions, the unsteady flow is governed by linear equations with variable 
coefficients which depend on the underlying steady flow. These equations are 
approximated using difference expressions determined from an implicit least-squares 
development and applicable on arbitrary grids. The resulting linear system of 
algebraic equations is block tridiagonal, which permits an efficient, direct (i.e., 
non-iterative) solution. 

In previous work unsteady solutions were obtained on a rectilinear-type, non- 
orthogonal, body-fitted, and periodic grid (termed "the cascade mesh") for cascades 
of sharp-edged airfoils aligned with the steady flow. Under the present effort the 
solution procedure has been extended to treat blades with rounded or blunt edges at 
incidence relative to the inlet flow. As part of this effort analytical studies 
have been conducted to clarify the behavior of first-order perturbation solutions 
in the vicinity of blade edges. Further, the numerical approximation has been 
applied to determine unsteady solutions for blunt leading-edge blades on the 
cascade mesh aind on a dense "local" mesh surrounding the leading edge. 

Numerical results are presented to Illustrate the influence of blade geometry 
(including leading edge curvature), inlet Mach number, and Inlet flow angle on 
unsteady response. Cascade and local mesh solutions indicate that leading edge 
curvature effects must be taken into account to accurately predict unsteady force 
and moment components in-phase with blade displacement. Predictions for staggered 
and unstaggered cascades of NACA 0012 and flat plate airfoils reveal that unsteady 
aerodynamic response is strongly dependent on blade thickness and inlet Mach 
number. In particular, present results indicate that the torsional stability 
margin for staggered cascades is substantially less for NACA 0012 than for flat 
plate blades. Further, the stability of out— of— phase bending or torsional motions 
is enhanced with increasing inlet Mach number. For subsonic attached flow, inlet 
flow angle or mean incidence variations appear to have only a minimal effect on 
unsteady force and moment. However, inlet flow angle variations do cause 
substantial changes in unsteady pressure distributions near blade leading edges. 
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INTRODUCTION 


A research program is being conducted to develop an aerodynamic analysis 
for predicting unsteady subsonic or transonic flow in a compressor or turbine 
blade passage. This analysis is intended to eventually serve as the aerodynamic 
component of a flutter or resonant stress design prediction system for turbo- 
machinery blade rows. To date, the effort has been concentrated on the subsonic 
flutter or free-vibration problem, but extensions of the analysis to treat 
transonic flows and forced aerodynamic excitations are planned as future work. 


Background 

Substantial progress has been achieved on the development of semi-analytical , 
subsonic and supersonic, unsteady cascade solutions (e.g. Refs. 1-5). However, 
these are essentially based on classical linear aerod 3 uiamic theory in which the 
unsteady flow is treated as a small disturbance about uniform steady flow. Such 
solutions do not apply if velocity gradients due to incidence, blade shape, or 
operation in the transonic Mach number regime are significant. Hence, they fail 
to meet the needs of turbomachinery designers over a wide range of practical 
operating conditions. 

To partially overcome the foregoing limitations an aerodynamic model has been 
formulated for subsonic inlet and exit conditions which includes the effects of 
blade geometry and flow turning on unsteady response (Refs. 6, 7). Here, the 
unsteady flow is considered as a small-amplitude, harmonic fluctuation about 
nonuniform steady flow. An asymptotic expansion of the time-dependent velocity 
potential then provides equations which govern the steady and small-disturbance 
unsteady flows in a single, extended, blade-passage region of the cascade. The 
steady flow is determined as a solution of the full-potential equation and the 
unsteady flow is governed by linear equations with variable coefficients which 
depend on the steady flow. Key features of this formulation which permit 
efficient numerical prediction of unsteady forces and moments are as follows: 
the unsteady boundary value problem is linear and time-independent; boundary 
conditions are applied at the mean position of blade and wake surfaces; the unsteady 
flow is periodic from blade to blade; and near-field numerical solutions are matched 
to far-field analytical solutions. Thus a numerical approximation is only required 
over a single blade-passage region of finite extent. 

Numerical solutions of the nonlinear steady problem are currently available 
for subsonic (Ref. 8) and transonic (Refs. 9, 10) flows. A numerical model for 
solving the unsteady equations on a rectilinear-type mesh which spans the extended 
blade passage region, but is not sufficiently detailed to permit an accurate 
resolution of the flow in the vicinity of blunt leading or trailing edges, has 
been developed as part of the present research program (Ref. 11). Unsteady 
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equations are solved on a nonorthogonal, body-fitted, and periodic mesh (herein 
called the cascade mesh) which facilitates the implementation of blade, wake 
and cascade-periodicity conditions, but disallows the use of standard difference 
approximations. Instead difference approximations are based on an implicit 
least-squares development applicable on arbitrary grids. This permits flexibility 
in the choice of difference neighbors and the simultaneous approximation of 
differential equation and boundary condition at boundary points, strategies which 
can significantly improve numerical accuracy. 

Unsteady response predictions based on the foregoing aerodynamic and numerical 
models have been reported (Refs. 7, 11) for subsonic cascades of sharp-edged 
blades with mean camber lines aligned with the steady flow. Response predictions 
for flat-plate cascades are in very good agreement with those based on Smith's 
analytical solution (Ref. 2) for both subresonant and superresonant blade 
motions (Ref. 7). Predictions for unstaggered cascades consisting of double- 
circular-arc or thin-circular-arc profiles reveal that blade thickness produces 
significant coupling between the steady and unsteady flows while steady blade 
loading (or flow turning) due solely to blade camber causes only weak steady/ 
unsteady interactions (Ref. 7). Recent results for staggered cascades of flat- 
plate and double-circular-arc blades (Ref. 11) Indicate that the coupling between 
steady and unsteady flows caused by blade thickness is particularly strong at 
high Mach number and at low vibration frequency. 

Similar linearizations with respect to a nonuniform steady flow have been 
considered in Refs. 12-16. Guiraud-Vallee, et. al. (Ref. 12) have studied 
compressible flow past an oscillating airfoil in a channel and have presented 
numerical results for subsonic and transonic flows past double-circular-arc and 
NACA 63A015 airfoils. Atassl and Akai (Refs. 13, 14) have derived analytic 
solutions for incompressible flow past a cascade of Joukowski airfoils. They 
have presented extensive results which indicate that blade thickness and incidence 
have a significant impact on unsteady response. Whitehead and Grant (Ref. 15) 
have developed finite-element steady and unsteady solution algorithms for subsonic 
cascade flows and have determined results for unsteady moments on a turbine 
cascade which are in agreement with experiment. Finally, Ni and Sisto (Ref. 16) 
regard the fluid properties as dependent variables and employ an explicit time- 
marching numerical procedure to determine solutions to the Euler equations of 
motion. This approach requires substantial computing time, but it is applicable 
to transonic cascade flows with strong shocks. 


Scope of the Present Effort 

It has been clear from previous work that improvements in the unsteady 
numerical approximation would be required before a practically useful flutter 
prediction scheme could be realized. Of most Immediate concern was the development 
of solution procedures for general blade profiles (i.e., rounded or blunt-edged 
blades) and the effects of mean flow incidence. Thus under the current phase of the 
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overall research program, a numerical approximation based on concepts from singular 
perturbation theory (Ref. 17) has been adopted for resolving the unsteady boundary 
value problem. Unsteady equations are first solved on the "rectilinear" cascade mesh 
to obtain an "outer" approximation to the linear unsteady flow, and then solutions are 
determined on dense local meshes to obtain "inner" approximations in limited regions 
of high velocity gradient. The solution on the cascade mesh provides outer boundary 
condition data for the subsequent local solutions. 

This procedure is applied here to determine the unsteady flow past cascades of 
blunt leading-edged airfoils. The unsteady numerical approximation on the cascade 
mesh has been modified so that errors caused by an inaccurate resolution of the flow 
near a blunt leading edge are confined to a limited region containing the leading edge. 
In particular, mesh points are not placed "too close" to the leading edge, and steady 
flow gradients are replaced by outer approximations in the vicinity of the leading 
edge. In addition, solutions on a dense local mesh which conforms to the shape of the 
airfoil and surrounds the leading edge are obtained using the implicit, least-squares, 
numerical approximation previously developed for the cascade mesh calculation. Unsteady 
response predictions for cascades of vibrating NACA 0012 and flat plate airfoils are 
presented to demonstrate the solution procedure and to investigate the effects of 
compressibility, blade geometry and mean incidence on unsteady response. 

As a part of the current effort an analytical study has been conducted for the 
limiting case of quasi-steady, imcompressible flow to gain a better understanding of 
the impact of leading and trailing edge geometry on the behavior of a first-order 
perturbation solution. Results of this study are described in detail in the APPENDIX 
to this report. They indicate, that for a blunt-edge airfoil, the asymptotic 
series used in the present aerodynamic model will converge to the exact solution 
for the flow, if the airfoil displacement is sufficiently small. Further, the first- 
order term of the asymptotic series will be analytic throughout the flow, and is thus 
susceptible to a numerical resolution, provided that a sufficiently dense mesh is 
employed in the vicinity of a blunt leading or trailing edge. 
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THE UNSTEADY AERODYNAMIC MODEL 


Basic equations which describe the flow past a finite-deflection cascade 
of airfoils undergoing small amplitude, harmonic vibrations are presented below. 
These equations apply to subsonic or transonic flows; i.e., flows with subsonic 
inlet and exit velocities, but with embedded supersonic regions adjacent to 
blade surfaces. The derivation essentially involves a perturbation expansion 
of the velocity potential which provides a nonlinear boundary value problem 
for the zeroth order or steady flow potential and a linear, variable-coefficient 
boundary value problem for the first order or unsteady potential. In the 
following discussion all quantities are dimensionless. Lengths have been 
scaled with respect to blade chord, time with respect to the ratio of blade 
chord to upstream free-stream speed, and pressure with respect to the upstream 
free-stream dynamic pressure. 


Description of the Problem 

Consider Isentropic and irrotational flow, with negligible body forces, 
of a perfect gas past a two-dimensional oscillating cascade (Fig, 1). The 
velocity potential, $(X, t), is governed by the time-dependent, full-potential 
equation; i.e., 

= 2V4> V4>,+ 74> V(7*) /2 (1) 

The speed of sound propagation, and the fluid pressure, are determined from 
Bernoulli's equation and the isentropic relations. It follows that 

= I-(7-I)M_1 {4>^ + [(7*f-|]/2} (2) 

where M is the Mach number, y is the specific heat ratio of the fluid and the 
subscript -“ refers to upstream free stream conditions. 

The blades of the cascade have finite thickness, camber and mean incidence 
angle relative to the inlet flow, and are undergoing identical, small amplitude, 
harmonic motions at frequency m and constant interblade phase angle, a. The 
mean or steady state positions of the blade chord lines coincide with the line 
segments mG^^ ^ ^ ^ 1 ^ ™ 0» ^ 1» i 2,.,., where and Gy are 

the X and y - components of the cascade gap vector, G, which is directed along 
the locus of blade leading edges with magnitude equal to the blade spacing 
(Fig. 1). Blade shape and orientation relative to the stream, and the amplitude 
and frequency of the blade motion are assumed to be such that the flow remains 
attached to the blade surfaces; i.e., 

(7*n)^ , 

m rrt 


m =0, ±l,± 2, . .. 


(3) 


Here(Rjj is a vector measuring the displacement of points on the instantaneous 
position of the mth blade surface, relative to their mean positions, t is 
time, and n is a unit outward normal vector. In addition to Eq. (3), the fluid 
pressure and normal velocity must be continuous across wake surfaces, and 
acoustic energy must radiate into the far field. In the inviscid approximation 
the wakes are represented by thin vortex sheets, each of which, emanates from 
a point in the vicinity of a blade trailing edge and extends downstream. 

To first order the relative displacement vector for the mth blade is given 
by 

(Rni=i'e =(h + axRp)e (A) 

where h defines the amplitude and direction of blade translations, a defines 
the amplitude of blade rotations, and Rp is a vector extending from the mean 
position of the mth axis of rotation (i.e., Xp + mG^, Yp + mG^) to points on the 
mean position of the mth blade surface, S^,. Only rigid motions are considered; 
however, the formulation could be readily extended to include elastic deformations 
of the blades. The components h^^, hy, and a of the vectors h and a are, in general, 
complex to permit phase differences between the translations in the x and y-direc- 
tions and the rotation. These rigid two-dimensional motions describe chordwise, 
bending, and torsional vibrations of actual rotor blades. 


Perturbation Procedure 

A perturbation approach serves to replace the nonlinear, time-dependent 
boundary value problem for the velocity potential, ^ (X, t), by two time-indepen- 
dent boundary value problems for a zeroth order or steady potential, (X), and 
a first order or unsteady potential, (J> (X) e^“^. In addition, unsteady boundary 
conditions can be applied at the mean positions of the blade and wake surfaces, 
and blade pressure distributions and aerodynamic response coefficients can be 
evaluated in terms of information prescribed at blade mean positions. The 
velocity potential is expanded in an asymptotic series in e; i.e., 

♦i(x’,t) = ♦ (x) + ^(x) e'"’ + Ck*) (5) 

j.O 

where X is the position vector and t is a small parameter related to the 
amplitude of the blade motion. In addition, Taylor series expansions; e.g., 

v4[ - v4| +(^^ v)v4| (6) 
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are applied to refer boundary conditions to the mean positions of the blade and 
wake surfaces. The mean positions, of the unsteady wakes are assumed to 

coincide with the steady flow stagnation streamlines. After substituting the 
foregoing series into the full governing equations, equating terms with like 
powers in e, and neglecting terms of second and higher order in e, nonlinear 
and linear variable-coefficient boundary value problems are obtained, respectively 
for the zeroth and first order potentials. 


As e 0 the blade and wake surfaces collapse to their mean positions. 
Hence, the zeroth order term in the power series expansion, Eq. (5), is the 
velocity potential, 4> (X), due to steady flow past a stationary cascade. The 
cascade geometry, the prescribed form of the blade motion, and the linearity 
of the first order equations, require that the first order or unsteady potential 
be harmonic in time (i.e., (X, t) = (f> (X) and that the steady and 

first order unsteady flows exhibit blade-to-blade periodicity; i.e. , 


♦ lx) =^(x + mG) - mV_„ G 
^(x) =^(x+mG)e“'*”‘^ 


m = 0,± I, ±2, ... 


(7) 


The periodicity conditions permit numerical solutions for the steady and unsteady 
flows to be determined in a single, extended , blade-passage region of the cascade 
(e.g., the region ^ in Fig. 2). Equations governing the steady potential, assumed 
to be known in the present study (Refs. 8-10), follow from the differential 
equation, Eq. (1), Bernoulli's equation, Eq. (2), and the flow tangency condition, 

Eq. (3), after replacing i' (X, t) by 4> (X) and omitting time derivative terms. 

Steady subsonic velocities are assumed to be approximately uniform beyond some finite 
distance from the blade row (e.g., or in Fig. 2); i.e., 




( 8 ) 


The Unsteady Boundary Value Problem 

It follows from Eqs. (1), (2), and (5) that the unsteady potential is governed 
by a linear differential equation of the form 








7 (^ 


(9) 


7 



(10) 


£? - i_ 
Dt ■ at 


+ V4> *7 


i ui + 7fl) *7 


and A is the local speed of sound in the steady flow. The influence of 
blade shape and steady turning appears in this equation through the variable 
coefficients (A^, (Y“1 )V^iJ>, etc.) which depend on spatial derivatives of the 
steady potential. Upon substituting the series expansions, Eqs . (5) and (6), 
into the flow tangency relation, Eq. (3), and noting the = ng+Ccixn*g) e^“*^ + 
O(e^), it follows after some algebra that the first order flow tangency condition 
has the form 


7(^-n = [iuir+a x7*-C7 •7)7<t>] * [itijr-7fr 74>)]- ne*'”®’ 

on S^, m = 0,± I, ± 2. . . . 

The first term on the right-hand-side of Eq. (11) is the surface velocity. 

The second term accounts for the effects of blade rotation relative to the 
mean flow and motion through a varying mean velocity field. The conditions 
of continuity of pressure and normal velocity across blade wakes and the blade- 
to-blade periodicity of the unsteady flow require that 

«o) 

) on W^,m = 0,±l,±2,... (12) 

A„{7^-n) =0 I 



where A (p denotes a difference across the roth wake; e.g., 
m 


An,<#> = 4> XonWm ( 13 ) 

The unsteady potential in the reference extended passage (i.e., the region 
bounded by the upper and lower surfaces of the zeroth and first blade, respectively, 
periodic lines upstream of the blade row, and by the zeroth and first steady flow 
stagnation streamlines downstream of the blade row (Fig. 2)) is determined as a 
solution of the variable coefficient differential equation, Eq. (9), subject to 
the flow tangency requirement, Eq. (11), at the upper, sj, and lower, S^, surfaces 
of the zeroth and first blades, and the continuity conditions, Eq. (13), across 
the zeroth wake, Wq. The unsteady velocity is continuous and differentiable 
upstream of the blade row, and hence, the blade-to-blade periodicity condition, 

Eq. (7) , is used in conjunction with the unsteady differential equation to 
determine values of the unsteady potential along upstream periodic boundaries. 
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Unstaady Far-Field Solutions 


To complete the specification of the unsteady boundary value problem conditions 
at upstream and downstream boundaries of the extended blade passage region must be 
established. In general, unsteady disturbances do not attenuate in the far-field 
of the cascade, and hence, it is difficult to place explicit conditions on the 
unsteady potential, cji , at upstream and downstream boundaries of the extended blade 
passage region. Instead, analytic far-field solutions, based on the uniform steady 
velocity approximation, Eq. (8), have been determined for subsonic inlet and exit 
conditions. These solutions can be matched to the near field numerical solutions 
at finite distances upstream and downstream of the blade row (e.g., on the lines 
5 = 5 in Fig. 2) . 

+ 

The unsteady potential is continuous far upstream of the blade row (i.e., 

4i = <t)j, for 5 < K_) and has both continuous and discontinuous components downstream 
of the blade row (i.e., = <1)^ + for 5 > 5+) • The continuous potential accounts for 

acoustic wave propagation into the far-field. Fourier methods can be used to provide 
expressions for the continuous components of the far-field potential in terms of 
the independent variables 5 and n where the 5, n-coordinate axes are perpendicular 
and parallel, respectively, to the cascade inlet plane (Fig. 2). It follows that 
(cf.. Refs. 6 and 7) 


j«»GO + + * 1 + 

where the coefficients, b. ,are obtained from the relation 

^ ' ^j.(C,Tj)exp (-iqj»>)dTj , 

The discontinuous component of the unsteady potential far downstream results from 
the counter vorticity shed from the trailing edges of the blades and convected 
along the blade wakes. A closed form solution (Refs. 6, 7) for this potential 
has been determined; i.e., 

F (u)exp (-twT/Vgj), €^€4- 


(14) 


(15) 


in terms of the Cartesian coordinates t and v, with the T-axis coinciding with the 
far-downstream, zeroth wake (Fig. 2). The constants q ,Xj _and the function 
F(u) in the foregoing equations are defined expllcity In Ref. 7. 
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Aerodynamic Response Coefficients 


Solutions to the steady and unsteady boundary value problems are required 
to determine blade pressure distributions and aerodynamic force and moment 
coefficients. The pressure at the mth blade surface is given by (c.f. Eq. (7)) 

^ P/ '^=0,±l,±2 (17) 


where P and pe^^^t ^re the zeroth (steady) and first (unsteady) order components 
of the fluid pressure, and S and S denote the moving and mean reference blade 
surfaces respectively. After expanding the velocity potential , $(X, t), in the 
manner indicated by Eqs. (5) and (6) and substituting the result along with 
Eq. (17) into Bernoulli's equation Eq . (2), it follows that 

9.-1/ ,2y/(y-i) 

P5 = 2(/M?„) '(M_„A)5 (18) 


and 






(19) 


Thus the steady and unsteady components of the pressure acting on a moving blade 
surface, ^ ^te evaluated in terms of information supplied at the mean position 
of the reference blade, S. The first term on the right-hand-side of Eq. (19) is 
the unsteady pressure at the mean position of the blade and the second is due to 
motion through a varying steady pressure field. The force and moment coefficients 
acting on the mth blade, i.e. 

A 

c. 




m = 0,± l,±2,. 


( 20 ) 


and 


C._= C. + Cm ('‘“'""'’•l (?<«'>. 


"m 


m= 0,± i,± 2,. 


( 21 ) 


are determined by simple integrations over the mean position of the reference 

blade. After some algebra it follows that the steady, Cp and C^, and unsteady 

c and c , force and moment coefficients are given by 
F M 

Cp=-^ Pg n ds, Pg Rp ds ( 22 ) 
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n ds + a X Cp , 




(23) 


where dS' is a differential vector tangent to the mean blade surface. It should 
be noted that the moment is taken about the moving pitching axis. The unsteady 
force and moment coefficients are the important results of an aerodynamic 
analysis intended for flutter predictions. Knowledge of these coefficients 
permits the evaluation of aerodynamic work per cycle and/or aerodynamic damping 
(Ref. 18) either of which can be used to determine whether the airstream supports 
or suppresses a prescribed blade motion. 
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THE UNSTEADY NUMERICAL MODEL 


In the approximate representation of a physical problem by a discrete 
numerical model, several components must be constructed. For the present case, 
these components include the discrete domain (or calculation mesh), the approxi- 
mating algebraic equations, and the solution procedure. Such components are 
generally not independent of each other with the choice of one influencing the 
choice of another. For instance, if the algebraic system is linear with proper 
structure, a direct solution procedure may be used, but if the system is nonlinear 
or linear but large and without particular structure an iterative technique may 
be necessary. It is important that the choice of one component not place overly 
restrictive constraints on the other components. 


Calculation Meshes 

Flow about a blunt edged airfoil has two length scales associated with it. 
The effects of airfoil geometry and motion on the undisturbed stream generally 
scale with blade chord and can be modeled on a calculation mesh of moderate 
density. However, flow phenomena near a blunt edge scale with edge radius and 
can be accurately modeled only on a much finer mesh. In this study, which 
involves flow through a cascade of airfoils with blunt leading edges and pointed 
trailing edges, the large scale effects of airfoil geometry and motion and their 
resulting flow are termed "cascade" while the small scale effects of the blunt 
leading edge and their resulting flow are termed "local". It would be difficult 
to approximate accurately and efficiently on a single mesh both the cascade and 
local effects because of the different length scales involved. Instead, in this 
study, the cascade flow is modeled on a cascade mesh of moderate density and the 
local flow is modeled on a local mesh of high density. 

The definition of an appropriate cascade mesh is much more difficult than 
the definition of an Isolated airfoil mesh. For the latter there are easily 
defined orthogonal meshes, e.g., potential-streamline meshes, "polar" meshes 
composed of normal and circumferential lines, and "parabolic" combinations. 

For a cascade these meshes are not so easily employed. For instance, the 
potential-streamline mesh is not periodic since potential lines generally are 
not parallel to the cascade. Boundary conditions arising from cascade periodicity 
in this case the periodic and wake boundary conditions of Eqs. (7,12) - can be 
approximated on such a mesh only with difficulty, since the mismatch of the 
mesh across the periodic boundary will require some sort of interpolation. 
Furthermore, the resulting linear system will not be tightly banded precluding 
the use of direct solution techniques. Orthogonal "polar" type meshes can be 
defined using conformal mapping as in Ref. 9. However, the necessary mappings 
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are not presently known for all cascades particularly those with high turnings 
or small gap to chord ratios. In contrast, the definition of an appropriate 
local mesh depends only on airfoil geometry and so is independent of whether the 
airfoil is isolated or placed in a cascade. 

The cascade mesh used in this study is the periodic and body fitted but 
nonorthogonal one shown schematically in Fig. 3. It is composed of "axial" 
lines, which are parallel to the blade row making the mesh periodic, and 
"tangential" curves, which are percentile averages of the upper and lower 
boundaries making the mesh body fitted. The leading and trailing edge axial 
lines are not in the cascade mesh eliminating the need for special considerations 
at leading and trailing edge points. The local mesh is the "polar" mesh shown 
schematically in Fig. 4. It is composed of "radial" lines normal to the airfoil 
surface and "circumferential" lines which roughly parallel the airfoil surface. 

Two Step Calculation Procedure 

The unsteady boundary value problem of Eqs. (9-16) is modeled on the 
cascade mesh using the algebraic approximations described below. In order to 
remove local phenomena Induced by a blunt leading edge from the cascade solution, 
steady potential derivatives are modified to remove leading edge curvature effects. 
This results in a cascade solution which is akin to an "outer" solution. These 
leading edge effects are then restored to the solution calculated on the local 
mesh with unmodified steady potential derivatives. Cascade effects are introduced 
to the local solution by interpolating the potential calculated on the cascade 
mesh to the outer radial and circumferential boundaries of the local mesh and 
imposing the interpolated potential distribution as a boundary condition. The 
solution to the unsteady boundary value problem is then taken to be the 
local solution in the local mesh region and the cascade solution elsewhere. 

To date, the cascade and local flows have been calculated by the two-step 
procedure described above, in which the local solution is essentially a correction 
to the cascade solution near the leading edge. Because the local leading edge 
effects are not in the cascade calculation, it is necessary to choose the local 
region extensive enough so that leading edge effects become insignificant on the 
local mesh boundaries. This insures that the local boundary value problem has 
valid information on the outer boundary. It is planned eventually to iterate 
between cascade and local solutions until the two solutions balance at the 
interface. When this is done a much smaller local region will suffice. 


Algebraic Approximations 

The nonorthogonality of the cascade mesh precludes the use of Cartesian difference 
approximations based on Taylor series expansions and so requires the development 
of more general approximations. These approximations are then used in both the 
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cascade and local problems. Consider a linear differential operator,^, which 
operates on a constant by multiplying that constant by qO. is to be approximated 
at the mesh point Qq in terms of the values of ij) at Qq and at certain neighboring 
mesh points, Q^, which, with Qq, are called a neighbor set. Quantities 

subscripted with 0 or m are evaluated at the corresponding neighbor set member. 

The approximation is developed in terms of an implicit interpolation of the form 



(24) 


where the are prescribed interpolating functions which vanish at Qq - homogeneous 
polynomials, for instance - and the are interpolating coefficients, which must 
be determined. The definition of y depends on the specific interpolation used, 
but, in any case, it is assumed to depend on values of (P at Qq and its neighbors; i.e. 


X 

where ^ =(5<})^, . . . , 5cj) ) , and the interpolating matrix, C , and vectors, and 
c^, depend on the geometry of the neighbor set but are independent of ip , In the 
present analysis the interpolation arises from a least squares fit, described 
below, and thus the method is termed an implicit least squares approximation. 
After replacing 61 }) by its interpolate. It follows that 


where is the Nxl vector with components ./f". Further, after combining 
Eqs. (25) and (26), the differential expression ^<p is approximated by the 
difference expression 


X + /8o^O + 

m" I 

where L is the difference operator derived from ^ and the coefficients of the 
difference operator are given by 

- I )9m+q0 (28) 

msi 
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3 being the Mxl vector with components 

in 

Since the differential operator Jt has no inhomogeneous term, (L(J>)q might be 
expected to have the form 


fn = l 

which emphasizes the "difference" aspect of its definition. Equations (27) and (29) 
would be equivalent if the vectors c® and c^ were zero. Although this is the case 
for most interpolations , c^ and c^ are nontrivial if the interpolate F is constrained, 

f\, '\i 

as it may be at a boundary where ip satisfies both a boundary condition (which is 
defined by the operator./ ) and the field equation (which becomes the constraint). 

Thus at a boundary, both the field equation and boundary condition can, in effect, be 
simultaneously approximated - a powerful feature of the implicit least squares 
approximation. 

Least Squares Difference Approximation 

The difference operator, L, which approximates the differential operator ./ 
is defined by Eqs. (27, 28) once the interpolation matrix and vectors are 
determined. First the interpolating functions, f’^, of Eq . (24) must be chosen. 

Common choices are the polynomials 

N=5« = f^=8T,. f^=8€St,, (30) 

2 2 

where , etc. Here, C and q could be any convenient 

coordinates but in this study, are taken to be the cascade axial and tangential 
coordinates Indicated in Fig. 2. The interpolation coefficients must then be 
defined. With eight neighbors but only five interpolating functions, it is 
impossible to define ^ such that F implicitly agrees with 6cJ) at all neighbors. 

On a rectangular mesh the reasonable but ad hoc Decisions that F agree with 6(|) 
at the neighbors above, below, right and left and that a certain combination of 
F values at the corner neighbors agree with the same combination of 6i)> values 
leads to an interpolation matrix which produces the familiar difference 
coefficients obtained by Taylor series expansions. On a general mesh these ad 
hoc decisions are no longer reasonable. What is needed is a formal method of 
defining the interpolating coefficients from the available data regardless of 
the geometry of the neighbor set. Least squares techniques provide such a 
method. 
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A measure of the error in replacing 6(() with its interpolate, F, in Eq. (24) is 


f 


8^] (31) 

where the over bar indicates complex conjugation and the Wj^ are chosen to give 
proper importance to the various neighbors. In practice w^ is usually given by 


Wm= >/|Qm-Qo* (32) 

but it may be set to zero if the corresponding neighbor is to have no influence. 

The real, nonnegative function e is minimized as a function of the real and imaginary 
parts of y. There results the following complex linear system to determine y : 

4y=^8^ (33) 

The NxN matrix /i and NxM matrix ^ are given in terms of the MxN matrix, ^ = (f^) = 
(f'^(Q^)) and the MxM diagonal matrix, W = diag (wj^, ..., w^) by 

( 34 ) 

where the superscript * indicates conjugate transposition. The interpolation 
matrix is then given by 


Q- (35) 

and the Interpolation vectors are zero. 

At a boundary point, Qq, the operator, J. , to be approximated usually is 
derived from a boundary condition rather than from the field equation, which 
probably also holds at Qq. Suppose the field equation, Eq. (9), is generalized 
to have the form 




( 36 ) 
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where J is a linear differential operator which operates on a constant by 
multiplying the constant by h®. Approximate A <t> by + F) and apply at 

Qq to obtain 


g(y) ■ = 0 


(37) 


where g(y) is a complex constraint functional and Af has components 


Equa 11011^(37) describes an algebraic constraint on the interpolate, F. The 
coefficient vector, 7 , is then determined by minimizing the error function, 
e, subject to this constraint. This is accomplished by using a complex Lagrange 
multiplier and minimizing the constrained error function 


ec(y A)= e(y) + Re{>} Re {g(y)}+lm |x}lm{g(y)} 


(38) 


with respect to the real and imaginary parts of 7 and the Lagrange multiplier , X . 
Equation (33) becomes 

Y +(o^f)Q X 




The interpolation matri:; and vectors are then given by 


where H is the scalar 


£o= - -V ^ 

Ci= ^ /3-'(>f)o 


(39) 


(40) 




(41) 


Let s be the scalar 




(42) 


Then, from Eq . (28) 




(43) 


Thus 6 is defined as if the operator J -s A viere being approximated; i.e., 

the operator being constrained is a particular linear combination of the boundary 

condition operator and the constraint operator. 
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Solution Procedure 


For the cascade calculation, the neighbor sets are generally defined as shown 
in Fig. 3; i.e., "centered" neighbor sets are defined if Qq is a field point, and 
"one-sided" neighbor sets are defined if Qq lies on a blade or wake boundary. 

In any case, all neighbors fall on the axial mesh line through Qq and on the two 
immediately adjacent axial lines. This placement is important to obtain a desirabl 
structure for the solution matrix. Derivatives of the unsteady potential 
exhibit singular behavior at sharp edges, in that near leading edges 

and near trailing edges where r^denotes distance from the edge (Ref. 19) 

Rapid, near-singular, variations are also observed near blunt edges. Thus for 
cascade calculations, either the first point(s) on the blade or the first point 
off the blade on the upper and lower boundaries of the passage (Fig. 2) Is deleted 
from neighbor sets near blade edges to ensure that differences are not taken 
across singularities. Affected neighbor sets are said to be adjusted. 

At field points the unsteady differential equation^ Eqs. (9,10), is 
approximated using centered neighbor sets. For points on the upstream periodic 
boundary (Fig. 2), values at neighbors above the mesh region are related to (|) 
values at points within the mesh region by the periodic condition, Eq . (7). On 
the lower periodic boundary the periodic condition is directly applied, since to 
approximate the field equation on both upper and lower periodic boundaries would 
be redundant. For points on the far-upstream boundary (C = C_) , <i> values at 
neighbors upstream of the mesh region are related to $ values on the boundary 
using the far field condition, Eq. (14), with the Fourier coefficients being 
evaluated using Eq . (15). The infinite sum in Eq. (.14) is truncated after a 
few terms - usually three- in each direction, and the integral in Eq . (15) is 
approximated using a trapezoid rule quadrature. A similar procedure is followed 
on the far-downstream boundary (C = , but here the discontinuous component of 

d given in Eq . (16) must be taken into account. At blade points the flow tangency 

condition, Eq . (10), is approximated using "one-sided" neighbor sets. At wake 
points the pressure continuity and normal velocity conditions of Eq. (12) are 
also approximated using "one-sided" neighbor sets. Letting Ly“ and represent 

either one of the approximating difference operators on the upper and lower 
wakes respectively, the resulting difference equations take the form 



where Qy+ is a mesh point on the lower wake and Q^- is the corresponding point 
on the upper wake. One of these (continuity of pressure) is applied on one 
wake boundary and the other (continuity of normal velocity) is applied on the 
other wake boundary. At both blade and wake points, the difference equations are 
constrained to satisfy the field equation. 
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For the local calculation matters are much simpler since no special 
considerations (e.g., neighbor set adjustment) need be taken at the leading edge 
and the many types of boundary conditions have become only two. Neighbor sets 
are chosen as for the global calculation, i.e., centered at field points and 
one-sided at boundary points. The unsteady differential equation, Eqs. (9, 10), 
and the flow tangency condition, Eq. (11), are approximated as above at field 
points and blade points respectively. At points on the outer boundary of the local 
mesh the unsteady potential is given the value interpolated from the cascade 
calculation. 

Let (() be a vector of the (p values on the £,th axial mesh line for the cascade 
calculation or on the fth radial mesh line for the local calculation. Because 
neighbor sets of points on the S,th line include points only from lines i-1, SL, 

SL+1, (and because the cascade calculation mesh is periodic,) the systems of 
linear algebraic equations for both cascade and local calculations have the 
following block tridiagonal form. 



+ C, <#>2 = d, 



= it . 2Si<L-l 

(45) 


--iL 



With this structure, the systems can be solved directly and efficiently using 
Gaussian elimination as described in Refs. 20 and 21. The sub-matrices A and 
C^, are sparse (as is B^^^) being basically scalar tridiagonal. This can be taken 
advantage of to improve computational efficiency during the forward elimination 
phase of the solution (Ref. 21). If a non-periodic mesh were employed for the 
cascade problem, the block tridiagonal structure would be lost and an iterative 
procedure would be required to obtain a solution. 
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NUMERICAL RESULTS 


The unsteady aerodynamic model accounts for the effects of blade geometry, 
including blunt leading and trailing edges, and mean incidence on unsteady 
response. However until recently, the unsteady numerical model was applied 
only on a rectilinear- type mesh which spanned the extended-blade-passage region 
(the cascade mesh shown in Fig. 3). Thus previous unsteady response predictions 
were restricted to cases in which the flow near blade edges could be reasonably 
approximated on this mesh; i.e., to cascades consisting of sharp-edged blades 
with mean camber lines aligned with the steady flow (Refs. 7 and 11). Predictions 
for flat plate cascades were shown to be in good agreement with Smith’s (Ref. 2) 
semi-analytical results for both subresonant and superresonant blade motions. 
Further predictions for double-circular-arc and thin-circular -arc blade profiles 
revealed that blade thickness produces a strong coupling between the steady 
and unsteady flows, particularly at high Mach number or low vibration frequency, 
while steady flow turning due to blade camber causes only weak steady/unsteady 
interactions . 

Under the present research program the numerical model has been developed 
further to treat cascades of blunt leading-edged airfoils at incidence relative 
to the inlet flow. The approximation on the cascade mesh has been modified so 
that errors due to an inadequate resolution of the large unsteady gradients near 
a blunt leading edge are confined to a limited region surrounding the leading edge. 
In addition, the numerical model has been applied on a dense local mesh (Fig. 4) 
to capture details of the rapidly varying flow in the neighborhood of a blunt 
leading edge. Thus for a cascade of blunt leading- edged airfoils, a numerical 
solution is first determined on a cascade mesh for the entire extended blade- 
passage domain, and then, a detailed solution is determined on a local mesh in 
the vicinity of the leading edge. The solution on the cascade mesh supplies 
boundary condition information (i.e., the unsteady potential distribution) on 
the outer boundary of the local mesh for the latter calculation. 


Present Examples 

The foregoing procedure has been applied to determine the unsteady flow 
past cascades of vibrating NACA 0012 airfoils. The mean position of the 
zeroth or reference blade surface is defined by the equation (Ref. 22) 


y(x)= ±5T[o.2969V)r-O.I260x-0.35l6x^ + 0.2843x3 - QI0l5x^],0< x < I 


(46) 
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where T = 0.12 is the blade thickness. The leading edge radius of the airfoil is 
1.1019 t 2 = 0.015867. For the present application the coefficient of the teirm 
has been changed to 0.1036 so that the upper and lower surfaces close at x = 1 
in a wedge-shaped trailing edge. Results have been determined for unstaggered 
(0= 0 deg.) and staggered (with 0= 45 deg.) cascades with unit gap-chord ratio 
(G = 1) for prescribed steady inlet Mach numbers, M_„, and flow angles, 
Corresponding exit Mach numbers. Moo, and flow angles, are determined by a 

global mass balance and the application of a Kutta condition, i.e., 

V ds I = ^•d’sj ^ , X = I (-^2) 

s s 


at blade trailing edges. Unsteady pressure difference distributions; i.e., 

Ap(x) = p-_(x)-p ^.(x)= p (x + G)e““^-p,+(x) (48) 

*» * , *0 


and aerodynamic response coefficients, Cj^ = cp.Sy and cj^, will be presented for 
cascades of NACA 0012 airfoils undergoing unit-frequency (o) = 1) , out-of-phase 
(a = * 180 deg.), single-degree-of-freedom, heaving (with hy = (1, 0)) or pitching 
(with a = (1, 0)) motions about an axis at mid-chord ((Xp, Yp) =(0.5, 0)). These 
rigid, two-dimensional, motions model bending or torsional vibrations of actual 
rotor blades. When the imaginary parts of the bending or torsional amplitudes 
are set equal to zero, the real and imaginary parts of the response coefficients, 
Ap, c^, and cj.j, are in-phase with the blade displacement and velocity, respec- 
tively. The stability of single-degree-of-freedom bending motions with Im {hy} = 

0 is governed by the sign of the Imaginary part of the lift coefficient. If Im 
(cl) < 0 the air stream tends to suppress the motion, and hence this motion is 
stable according to linear theory. Similarly, if Im {a} = 0, single-degree-of- 
freedom torsional motions are stable when Im {cj^} < 0 (Ref. 23). In addition 
to the NACA 0012 cases, for purposes of comparison, numerical results will be 
presented for oscillating flat plate cascades with blade mean positions aligned 
parallel to the free-stream direction. In this case $ = x and the unsteady 
equations reduce to those of classical small-disturbance theory (Ref. 19). 

Unsteady solutions for NACA 0012 and flat plate cascades were determined 
for a blade passage region extending two axial chords upstream and downstream 
of the blade row with far-field Fourier expansions, Eq. (14), being truncated 
after j = 3. For the NACA 0012 cases unsteady solutions were calculated on 
both cascade (Fig. 3) and local (Fig. 4) meshes. In addition, steady flows 
have been calculated on similar but coarser meshes using the finite area 
approximation developed by Caspar, Hobbs, and Davis (Ref. 8), and steady 
data has been defined on the unsteady meshes using constrained, polynomial, 
least-squares interpolations in an explicit manner. For the flat plate 
cases the steady flow is known and unsteady calculations are only required on 
the cascade mesh. 
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The unsteady results reported here were determined on cascade meshes con- 
sisting of 15 tangential lines and from 69 to 73 axial lines, 25 to 28 of which 
intersect blade surfaces. To minimize the effects of local errors on the cas- 
cade mesh solution, axial mesh lines did not Intersect blade leading and trailing 
edges. Further, steady derivative behavior in the vicinity of a blunt leading edge 
was modified to depend on a position vector emanating from the leading edge 
rather than on a position vector emanating from the center of the leading edge 
circle (c.f. the APPENDIX to this report) using an extrapolation based on the 
Prandtl-Glauert transformation. The spacing of the axial mesh lines as a 
percentage of axial chord varied from 20 percent at far-field boundaries to 
approximately 2 percent near blade edges to 5 percent near midchord, while the 
spacing between the tangential mesh lines as a percentage of distance between 
the upper and lower, extended-passage, boundaries varied from 2 percent near 
the boundaries to 15 percent near the middle of the passage. 

The local mesh employed for the NACA 0012 examples consisted of 12 circum- 
ferential lines which wrapped around the leading edge of the airfoil and 
extended from normal lines to the upper and lower surfaces at midchord, and 
100 "radial" or normal lines which extended outward from the airfoil to one- 
half the minimum distance (i.e. the throat) between adjacent blades. The 
normal lines were spaced so that circumferential distances half way out (i.e., 
at 25 percent of the throat) were equal. The distance between circumferential 
lines as a percentage of the throat varied from 0.125 percent near the airfoil 
to 10.5 percent at the outer boundary of the local region. Once steady derivative 
information has been stored, unsteady solutions on the cascade and local meshes 
described above require approximately 30 and 20 seconds, respectively, of CPU 
time on a UNIVAC 1100/81 system. It appears that the cascade and local meshes, 
used presently, are more dense than necessary for an accurate resolution of the 
unsteady flow, and hence, substantial reductions in computing time should be 
realized in future studies. 


Unstaggered Cascades 

Results for unstaggered NACA 0012 and flat plate cascades are shoxm in Figs. 
5 through 13. For the NACA 0012 cascade, inlet Mach numbers, M_oo,of 0.3, 0.5, 
and 0.6 and inlet flow angles, n_„, (Fig. 1) of 0 deg., 5 deg., and 10 deg. have 
been considered. For the flat plate cascade the inlet Mach numbers are 0.3, 0.5, 
0.6, and 0.7 and the inlet flow angle is 0 deg. in each case. Cascade and 
local mesh solutions for unsteady pressure difference distributions near the 
leading edge of the reference NACA 0012 blade are shown for out-of-phase bending 
and torsional vibrations in Figs. 5 and 6, respectively, for M_„, = 0.5 and 
n_oo = 0 deg. Here the pressure difference is plotted versus the square-root 
of the distance along the airfoil chord from the leading edge. Solutions on 
the cascade mesh approximate analytical singularities in Ap at the leading edge. 
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which are removed by the subsequent local solution. The latter determines a 
zero pressure difference at a blunt leading edge and local extrema In unsteady 
pressure difference near the leading edge. For the bending motion (Fig. 5) there 
Is a local minimum In Re {Ap} at x 0.0035 and In Im {Ap} at x 0.013, while 
for the torsional motion (Fig. 6) the real and imaginary pressure differences 
reach local maxima at x 0.0009 and x 0.012, respectively, and the real 
component has a local minimum at x 0.03. It should be noted that the singular 
behavior predicted by the cascade mesh solution for the real component of the 
pressure difference due to torsion (Re {Ap} -*-<*> as x ->■ O"*”; c.f. Fig. 6) differs 
from that predicted by classical linear theory (Re {Ap} as x O'*"; c.f. 

Fig. 11). This difference is due to the contribution of higher order singular 
terms (i.e., (Re {p} O (x~^^^) as x -»• 0) which appear in an "outer" approx- 
imation to the unsteady flow near a blunt leading edge. As indicated in Figs. 

5 and 6, significant differences exist between the real components of the lift 
and moment coefficients predicted by the cascade mesh solution and the combined 
(i.e., cascade plus local) solution. In the latter case local mesh results 
apply in the region covered by the local mesh and cascade mesh results apply 
elsewhere. However, the cascade solution provides a good approximation to the 
imaginary components of the response coefficients. 

The effect of inlet Mach nunber variation on the unsteady response to a 
cascade of vibrating NACA 0012 airfoils is illustrated in Figs. 7 through 9 for 
inlet Mach numbers of 0.3, 0.5, and 0.6. In each case the inlet flow angle is 
0 deg. Steady Mach number distributions on the reference blade surface are 
shown in Fig. 7. For the prescribed inlet conditions surface Mach numbers reach 
maximum values of 0.379, 0.667, and 0.864 at x = 0.132, x = 0.152 and x = 0.195, 
respectively. Results depicted in Figs. 8 and 9 indicate that the unsteady 
response is strongly dependent on inlet Mach number, particularly at high Mach 
number. For out-of-phase bending motions (Fig. 8) the imaginary component of 
the unsteady pressure difference increases near the leading edge, but decreases 
over the remainder of the airfoil as inlet Mach number increases. This behavior 
results in a decrease in imaginary lift component; i.e., an increase in aero- 
dynamic force opposing the blade motion. For out-of-phase torsional vibrations 
(Fig. 9), as M_oo increases the imaginary component of the unsteady pressure difference 
increases substantially over the forward 40 percent of the airfoil and decreases over 
the remaining 60 percent resulting in an increase in the unsteady moment opposing 
the motion (i.e., Im {cm> decreases). Hence, the effect of higher inlet Mach number 
is stabilizing for the two motions. 

Similar results are shown in Figs. 10 and 11 for flat plate cascades. Here 
the steady Mach number is a constant. An increase in Mach number has a stabil- 
izing effect on the motions of a flat plate cascade, but Mach number variation 
has a less pronounced effect on the aerodynamic response to vibrating flat 
plate blades than it does for the NACA 0012 profiles. A comparison between 
unsteady predictions for the NACA 0012 and flat plate cascades for a given inlet 
Mach number reveals the important influence of blade geometry on unsteady aero- 
dynamic response. The results provided in Figs. 8 through 11 indicate that the 
effect of blade thickness is stabilizing for the subsonic attached flows con- 
sidered, particularly at high Mach number. 
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The effect of mean incidence on the response to single-degree-of-freedom 
motions of unstaggered NACA 0012 cascades has been examined by computing 
unsteady solutions for an inlet Mach number of 0.5 and inlet flow angles of 
0 deg., 5 deg., and 10 deg. The numerical results indicate that an inlet flow 
angle variation from 0 to 10 deg. has little impact on overall unsteady pressure 
difference distribution, and hence, unsteady lift and moment coefficients. An 
Increase in mean incidence angle caused significant variations in unsteady sur- 
face pressure near the leading edge of the airfoil, but produced only a slight 
increase in stability margin for the bending or torsional vibrations. Real and 
imaginary unsteady surface pressure distributions near the leading edge of the 
reference blade are shown in Figs. 12 and 13, respectively, for the torsional 
case. Although there are substantial differences in unsteady pressure in the 
immediate vicinity of the leading edge (0 ^ Vx ^ 0.2) these do not persist over 
the entire airfoil, and therefore unsteady force and moment coefficients are 
relatively unaffected by the changes in inlet flow angle. 


Staggered Cascades 

Results for staggered cascades of NACA 0012 and flat plate airfoils are 
shown in Figs. 14 through 18 for a stagger angle of 45 deg. and inlet Mach numbers 
equal to 0.3, 0.5 and 0.7. The inlet flow angle is 5 deg. for the NACA 0012 
and 0 deg. for the flat plate cases. Surface Mach number distributions for 
NACA 0012 profile are shown in Fig. 14. Peak Mach numbers on the suction sur- 
face are 0.379, 0.649, and 0.970 and occur at x = 0.053, x = 0.068, and x = 0.095, 
respectively. The steady flows stagnate on the lower surface of the airfoil 
within a distance of 0.03 percent of blade chord from the leading edge and the 
stagnation point moves closer to the leading edge as inlet Mach number increases. 
Unsteady response predictions for out-of-phase bending (Fig. 15) and torsional 
(Fig. 16) vibrations of the NACA 0012 profiles follow similar trends with Mach 
number, M_co, as do the corresponding results for unstaggered cascades (Figs. 

7 and 8) . For the bending case the real component of the pressure difference 
distribution decreases substantially over the forward half of the airfoil leading 
to a reduction in real lift component and an increase in the real counterclockwise 
moment. The imaginary component of the pressure difference distribution increases 
near the leading edge, but decreases over most of the airfoil providing a sub- 
stantial increase in the magnitude of the force opposing the motion and a 
reduction in counterclockwise moment at high Mach number. For the torsional case 
real pressure differences decrease over most of the blade providing a reduction 
in lift and an Increase in counterclockwise moment. Further, imaginary pressure 
differences increase on the forward 30 percent of the blade and decrease over 
the remainder of the blade, leading to a reduction in lift and an increase in the 
moment opposing the motion. Hence, the effect of increasing inlet Mach number 
is to enhance the stability of both motions. 
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Results for the staggered flat plate cascades reveal a reduction in the 
imaginary component of the unsteady pressure difference due to bending (Fig. 

17) over most of the airfoil, and hence, an increase in the unsteady force 
opposing the motion with increasing Mach number. Further, there is a substantial 
increase in Im {Ap} due to torsion (Fig. 18) over the front-half of the airfoil, 
and hence, an increase in the moment opposing the motion, as Mach number increases. 
A comparison of the results given in Figs. 15 through 18 indicates that for the 
subsonic Mach number range considered, blade thickness enhances the aerodynamic 
stability of the bending motions, but causes a reduction in stability margin 
for the torsional motions. 

Unsteady Moment Due to Torsion 

In an effort to briefly examine the effect of blade geometry on overall 
torsional stability, limited additional calculations have been performed for the 
staggered (G =1, 0= 45 deg.) cascades of NACA 0012 and flat plate airfoils 
undergoing single-degree-of-freedom torsional vibrations about midchord. For 
these cases the inlet Mach number was set at 0.7 and the inlet flow angle was set 
at 5 deg. for the NACA 0012 cascade and at 0 deg. for the flat plate configuration. 
Results were first determined for unit frequency motions over the entire interblade 
phase angle range (i.e., -180 deg. < a < 180 deg.), and then for motions at 
to= 0.5 over the subresonant range of positive Interblade phase angles. At unit 
frequency, superresonant conditions (i.e. acoustic waves persist far upstream 
and/or far dovmstream of the cascade, c.f. Ref. 7) exist for -28.6 deg. 

< a < 112.4 deg. for the NACA 0012 cascade and for -29.4 deg. < a < 107.3 deg. 
for the flat plate cascade. The results for the imaginary component of the 
moment coefficient are depicted in Fig. 19. For superresonant motions the 
aerodynamic damping (which is proportional to -Im|cjv^f) is greater for the NACA 0012 
than for the flat plate blades. However, this situation is dramatically reversed 
over most of the subresonant region. The maximum calculated values of Im | Cj.j | 
were -0.157 at a = 120 deg. (subresonant) for the NACA 0012 cascade and -0.408 
at o = 60 deg. (superresonant) for the "classical" flat plate configuration. 

Hence, at least for this example, the effect of blade geometry is to cause a 
substantial reduction in torsional stability margin. In addition, the predic- 
tions at the lower frequency (co=0.5) indicate that the NACA 0012 cascade will 
experience subresonant torsional instabilities (according to linear aerodynamic 
theory) at low vibration frequencies. 


25 



CONCLUDING REMARKS 


The aerodynamic analysis described here has been developed to determine the 
flow past a finite-deflection, oscillating cascade operating in subsonic free- 
stream conditions. This analysis is intended to serve as a basis for investigating 
flutter or free vibration phenomena in jet engine fan, compressor and turbine blade 
rows. The unsteady aerodynamic model assumes small-amplitude fluctuations about a 
nonuniform steady flow and fully accounts for the effects of blade geometry and 
steady flow turning on the response to an oscillating cascade in fully-attached 
subsonic flow. The steady flow is governed by the full-potential, boundary value 
problem and the unsteady flow is governed by a linear boundary value problem with 
variable coefficients defined by the steady flow. With the present formulation 
a numerical resolution of time-independent , linear, unsteady equations is required 
only over a single extended blade-passage region of finite extent. Further, an 
incompressible model problem study (cf. the APPENDIX) indicates that the first- 
order or unsteady flow is analytic, but rapidly varying, near a blunt leading or 
trailing edge, and is thus susceptible to numerical resolution provided that a 
sufficiently detailed computation grid is employed in such regions. 

The unsteady numerical model is general and comprehensive enough to approximate 
accurately the unsteady equations despite the difficulties posed by cascade 
periodicity, nonrectangular geometry, and limited regions of rapid flow variation. 
Differential quantities are approximated in the physical plane using finite 
differences determined by an Implicit least-squares interpolation and applicable 
on arbitrary grids. The approximation is flexible enough so that neighbor sets 
can be adjusted to accommodate singular behavior near blade edges, and both 
differential equation and boundary condition can be simultaneously approximated 
at boundary points. In general, the unsteady equations are first solved on a 
nonorthogonal , body-fitted, and periodic cascade mesh (Fig. 3), and then on dense 
local meshes in regions of large velocity gradient (e.g., near a blunt leading 
edge. Fig. 4). On either cascade or local mesh, the unsteady boundary value problem 
is approximated by a block-trldiagonal system of linear algebraic equations which 
is solved by direct matrix inversion . 

In previous studies unsteady solutions were determined on a cascade mesh for 
cascades of sharp-edged blades with mean camber lines aligned with the steady 
flow (Refs. 7, 11), Predictions for flat plate cascades have been shown to be in 
good agreement with Smith's (Ref. 2) analytical results for both subresonant and 
superresonant blade motions. Further, predictions for double-circular-arc blades 
revealed that blade thickness produces a strong coupling between steady and 
unsteady flows, particularly at high Mach number or low vibration frequency. 

Under the present effort the solution procedure has been extended and applied 
to cascades of blunt leading-edged airfoils at Incidence relative to the inlet 
flow. Results have been presented here for cascades of NACA 0012 airfoils to 


26 


demonstrate solution capability and to further illustrate the influence of blade 
geometry (including leading edge curvature), Mach number, and incidence on unsteady 
response. Cascade and local mesh solutions indicate that leading edge curvature 
effects must be taken into account to accurately predict the unsteady force and 
moment components in-phase with blade displacement. Detailed local solutions for 
the NACA 0012 profiles provide unsteady pressure difference distributions which 
are zero at the leading edge and reach maximum or minimum values at distances 
from the leading edge on the order of the leading edge radius. Results for staggered 
and unstaggered cascades of flat plate and NACA 0012 airfoils reveal that unsteady 
aerodynamic response is strongly dependent on blade thickness distribution and 
inlet Mach number. In particular, the stability margin for torsional motions of 
staggered ( 0 = 45 deg. ) cascades is substantially less for the NACA 0012 than 
for the flat plate cases. The stability of out-of-phase bending or torsional 
vibrations is enhanced with increasing Mach number for both NACA 0012 and flat plate 
cascades, but particularly for the former. Finally, for subsonic attached flow 
past NACA 0012 cascades, mean incidence appears to have only a slight effect on 
unsteady force or moment. However, inlet flow angle variations do cause substantial 
changes in unsteady pressure distribution near blade leading edges. 


Future Studies 

Further developments in the unsteady aerodynamic analysis will be required 
before a comprehensive flutter prediction scheme can be realized. For the sub- 
sonic Mach number range in which blade vibrations are of practical concern, it 
will be necessary to extend the present analysis to include transonic flows with 
weak shocks. For transonic applications, shock motion conditions (Ref. 24) must 
be incorporated into the aerodynamic model and appropriate differencing procedures 
for supersonic regions must be determined for the unsteady numerical model. In 
addition, for vibrating thin fan or compressor blades at incidence - the positive 
incidence flutter problem (Ref. 25) - it appears that the aerodynamic model 
will have to be modified to incorporate the effects of viscous-induced, leading- 
edge, separation bubbles. Progress in this direction is heavily dependent on 
advances in the calculation of steady cascade flows. Finally, the aerodynamic 
analysis developed to predict cascade flutter or free vibration phenomena should 
be combined with solutions of the vorticity and entropy transport equations 
(cf. Refs. 26, 27) to provide an analysis for the prediction of forced excitation 
phenomena in turbomachinery blade rows. 
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APPENDIX 


INCOMPRESSIBLE MODEL PROBLEM 


Certain restrictions are inherent to a perturbation approximation. If these 
are not satisfied throughout the solution domain, the perturbation expansion is said 
to be singular, i.e., the asymptotic series does not provide a uniformly valid 
approximation to the true solution (Ref. 17). In classical linear theory this sit- 
uation is manifested by analytical singularities (i.e., |V(j)|''/ where r^. is 

the distance from the edge) in the first order solution at airfoil leading (n=-l/2) 
and trailing (n = 1/2) edges (Ref. 19). To clarify both the convergence properties 
of the asymptotic series associated with the present aerodynamic model (Eq. (5)) 
and the behavior of the first-order solution, especially in the vicinity of blunt 
leading edges, a quasi-steady (u^) , incompressible (M_co-K)) flow has been studied. 

The results of this study provide important considerations for a successful numerical 
resolution of the first-order unsteady flow field. 

Consider incompressible flow past an elliptic airfoil, M , which has undergone 
a rigid displacement from an original or mean position, S (Fig. 20). The latter 
is situated at angle of attack A relative to the uniform stream, V^. Body-fixed « , 
Y - and x,y-coordlnate axes coincide with the major and minor axes of the ellipses 
^ and S, respectively. Hence, ^ = (z-h)e~^“ where ^ =«+ , z = x + iy, h = h^ + 

ihy defines the amplitude and direction of the translational displacement, and a 
defines the amplitude and direction of the rotational displacement about ^ = 0. The 
complex potential, W, for the flow about the displaced ellipse, ^ , is given by 
(Refs. 28, 29) 

(l-b) W(^) = h{^,a) 

- F(^)cos A - i G(^) sin a + i (l-b) Vim, 




40 ‘ 


(Al) 


where 


^ ,G{^) = -v^^-4a^‘ 


r = -(i + b)sin A/2 


(A2) 


A 2 , 1/2 -' 

A = A + a, b is the thickness of the airfoil, a=(l-b) /4, and F is the circu- 

lation. The circulation has been specified by requiring that the rear stagnation 
point occurs at the trailing edge (^ = 1/2). For b = 0, Eq. (Al) defines the complex 
potential for the flow past a flat plate at angle of attack, A , relative to the 
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uniform stream V„. In this case the prescribed circulation, Eq. (A2) , guarantees 
that the flow does not stagnate near the trailing edge of the airfoil. 


Perturbation Solution 

The complex potential can be formally expressed as a power series in the param- 
eter e; i.e. , 


A OD i 

w(^) =2; 

j=0 


Wj(Z) = 


W{z)+ w{z) + (3(e2) 


(A3) 


where e depends on the displacement of the airfoil S from its original or mean posi- 
tion. The first two terms of this series are determined explicitly by expanding 
the second term on the right-hand-side of the expression 




(z^- 4a^) 


z^- 4a^ 


(AA) 


in a binomial series; i.e.. 


I + 





2(z2-4q2) 




(A5) 


substituting this result along with the series expansions for cosa and sina into 
Eq . (Al) , and neglecting second and higher order terms in hj^, hy , and a. After 
some algebra it follows that 


(l-b) w(z) = h(z, a) 


(A6) 


and (l-b)w(z) = a - s 


COS A 

(A7) 

+ s j^G(z)cos A - i{F(z)-(i-b^)/ 2 )sin a] /^z^- 4 q^ 

where s = h + iaz. Thus, W(z) is the complex potential of the flow past the ellipse 
S and w(z) is the linear (in a, h^^ and hy) approximation to the change in this solu- 
tion due to a rigid-body displacement. 


= - a F (z)sin A - ia 


G(Z) + 


(l-b^) 


z + 




4 q‘ 
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The power series, Eq. (A3) evaluated at the point z, converges to the exact 
solution, W, evaluated at ^ = (z-h)e~^“ for flow past the displaced ellipse, S , if 


= |(z-h)^ e |z*-40^| = |z + 2q|| 2- 2o| (A8) 

at each point z on or exterior to S. The foregoing inequality can be expressed in 
the form 


' ' 2 Z 


(A9) 


Thus, the power series, Eq . (A3), converges if the distance from z to ^ is small 
compared to the ratio of the product of the distances from z to the foci (at z = ± 
2a) of the ellipse S to twice the distance from z to the origin (z = 0) . In par- 
ticular, at the edges of the airfoil (z = ± 1/2), the foregoing inequality reduces 
to 


h ± ia/a] < b^4 - (?(lsl^) 


(AlO) 


2 

where b /4 is the edge radius. Thus, for sufficiently small displacements, the power 
series, Eq . (A3), will converge throughout the flow field, and the first two terms 
of this series will provide a valid first-order approximation to the flow past the 
displaced ellipse. If condition (AlO) is not satisfied (e.g., for a flat plate air- 
foil), the perturbation approximation is singular near the edges of the airfoil. 

The zeroth-and first-order complex velocities are given by 


= -^=-(z*-4o*) ^ [g{z) cos A -i [F(Z)-(l-b*)/z] sinA (All) 


-'0 = ias''(z-4aV'" [g(z)a F(z)- O-bh/z] 

-s(z^-4a^) [f(z)cos a -i G(z) sin a] 

-sz(z^-4a*) [g(z) cos a - i[F(z) - (i-b*)/ 2 ]sin aJ 


(A12) 


30 


For flow past a displaced ellipse the zeroth-and first-order complex velocities 

^ —1 /2 — /2 
are analytic. They behave like multiples of Zq ' and Zq ' > respectively, near the 

blunt leading and trailing edges of the airfoil, where Zq is a complex coordinate 

measured from the center of the edge circle. For a flat plate airfoil (b=0, a =1/4^ 

the zeroth order complex velocity has a square-root singularity (dW/dz<\j ^?(z+l/2)“^^^) 

at the leading edge and behaves like a multiple of (z - 1/2) at the trailing 

edge. In addition, the first-order velocity is singular like a multiple of 

(z + 1/2)”^'^^ at the leading edge and like a multiple of (z - 1/2)“^' at the 

trailing edge. For a flat plate at zero mean angle of attack (the classical linear 

problem) dW/dz = 1, and the order of the singularities in the first-order complex 

velocity is reduced by one (i.e., dw/dz (z + l/2)~^^^ at the leading edge and 

dw/dz (z - at the trailing edge). 


Governing Boundary Value Problems 

Since the zeroth order term, W(z) of the series expansion, Eq . (3), is the com- 
plex potential for the flow past the ellipse S (Fig. 1), <t> = Re {W(z) } is a solution 
of Laplace's equation which satisfies the following boundary conditions 


V (J> ■ "n = 0, on S 


(A13) 


and 


lim V «t> = cos A 6x + sin A 6y (A14) 

|zhco 

where n is a unit outward normal vector and e^ and 6y are unit vectors in the 
coordinate directions. The first-order velocity potential, (j) = Re {w(z)}, is also a 
solution of Laplace's equation, and it can be verified that ()>(x,y) meets the following 
conditions 


• n = X V<ti - (T-V) V4>j • n, on S 


(A15) 


and 


lim = O (A16) 

|zl-oo 
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f 


where r = (h^ + ay) e^ + (h - ax) ey. Thus exact solutions to the steady and the 
linear unsteady boundary value problems described in this report (c.f. Eqs. (1), 
(3), (8), (9) and (11)) have been determined for the flow past an isolated elliptic 
airfoil in the limiting case 0 and ui -»■ 0. Note that the unsteady boundary 

conditions of continuous normal velocity and pressure across the wake (Eq. (12)) 
are automatically satisfied since Vij> is continuous throughout the flow. Further, 
the unsteady far -field conditions (Eqs. (14) and (16)) are replaced by the condi- 
tion that the first-order velocity vanishes far from the airfoil. The zeroth-and 
first-order pressures acting on the displaced airfoil, S , follow from Eqs. (2), 
(18) and (19) for ^ 0 and w 0 and are given by 


Ps= i-(V<l>)s 


(A17) 


and 


Pj = [-2V4) - V<^ + 7 - VP]^= [-2V<t • V(^-7 V(V<t>)^]^ (A18) 


The zeroth-and first-order solutions, Eqs. (A6) and (A7),have been programmed 
and sample results for the first-order potential and normal velocity on the mean 
surface, S, and the first-order pressure difference are shown for rotational dis- 
placements with a = 2 in Figs. 21 through 25. Surface velocity potential distribu- 
tions for flat plate (b = 0) and elliptic airfoils (b = 0.1 and b = 0.2) at zero mean 
angle of attack. A, are depicted in Fig. 21. In each case the potential varies 
gradually over most of the airfoil, but abrupt changes in potential occur near the 
leading and trailing edges of the elliptic airfoils. The magnitude of these rapid 
changes increase with decreasing thickness (for b 0) and with increasing magni- 
tude of angle of attack. In contrast, abrupt changes in surface potential do not 
occur for a flat plate airfoil. 

Surface normal velocity distributions, as determined from the first-order flow 
tangency condition, are shown in Figs, 22 and 23 for flat plate and elliptic air- 
foils at 0 deg. and -10 deg. mean angles of attack. Details in the vicinity of the 
leading edge are shown in the latter figure. Results for the flat plate at A = 0 
(i.e., V(j) • ^Ig = ± 2) correspond to the classical-linear-theory boundary condition. 
The effects of thickness or non-zero mean angle of attack is to cause large and 
rapid variations in first-order normal velocity in the vicinity of airfoil edges. . 
These variations increase with decreasing thickness for b 0 or increasing magni- 
tude of angle of attack. For the flat plate at A = -10 deg. the surface normal 
velocity distribution is singular at the leading edge and varies monotonically with 
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distance from the leading edge along the upper and lov'er surfaces of the airfoil. 
However, for the elliptic airfoil the variations in first-order normal velocity 
along the upper and lower surfaces change sign near the edges of the airfoil (Fig. 

22) . This behavior has an Important impact on the numerical resolution of the first- 
order flow and will be discussed in more detail below. 

First-order pressure-difference distributions, Ap(x) , due to rotational dis- 
placements with a = 2, for symmetric mean flows (A= 0) past flat plate airfoil 
and elliptic airfoils with b = 0.1 and b = 0.2, and for flow past an ellipse with 
b = 0.1 at -10 deg. mean angle of attack are shown in Fig. 24. Leading edge details 
are provided in Fig. 25. The results depicted in Fig. 24 suggest that the first- 
order loading (i.e., Ap or Cj^) is more dependent on airfoil thickness than on mean 
angle of attack. Significant variations in the first-order pressure difference 
occur near the edges of the airfoil, particularly near the leading edge. For the 
blunt-edged airfoils the pressure difference distributions reach minimum values at 
a distance of b^/4 from the leading edge; i.e., at a distance equal to the leading 
edge radius (Fig. 25). 


Inner and Outer Expansions 

The results depicted in Figs. (21) through (25) indicate that flow parameters 
vary sharply in the vicinity of airfoil edges and gradually over the remainder of 
the chord. This behavior suggests that the flow past an elliptic airfoil depends 
on two length scales — the edge radius, b^/4, and the square root of the product of 
the distances from a given point, z, to the leading and trailing edges; i.e., 

(z^ - 1/4)^^^. Approximations to the zeroth-and first-order complex velocities, 
based on this concept, can be determined. An "inner" solution valid in a neighbor- 
hood of the leading edge, \ zq \ (? (b^) where Zg = z + 2a, is obtained by substi- 
tuting the expansion 


(z2-4a^)"~ (-yi-b2 )" [i-nzcVi-b^ (ai9) 


into Eqs. (A2) , (All) and (A12) . It follows after some algebra that for 
A -t (?(b)<<! 


(i-b) 


dw 

dz 


I 


(A+ib/2)Zc‘''^ + <?(b*) 


(A20) 


and 


O-b)^ 


s( A + ib/2) z^* /2 - a z 


- 1/2 




(A21) 
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Thus, as noted previously, the zeroth-and first-order complex velocities exhibit a 
and a behavior, respectively, in the vicinity of a blunt leading edge. 

Similarly, an "outer" approximation, valid in the region | z^ - 1/4 1 '\/^(l), 
is determined by substituting the expansion 


(Z^-4Q^)" ~ (Z^- 1/4)” ^l + nb*(Z*- 1/4) /4j + (b^) {klT) 

into the exact expressions for the zeroth-and first-order complex velocities. It 
follows for A (?(b)<<! that 




i + (z^-i/4)'^^ [iA(z-i/2)-bz] +^(b^) 


(A23) 


and 


(i-b) 


dw 

dz 


0 


-a[A + i(i+b)][i-U- 1/2)*^^ {z+i/2)'^^] 


(A24) 


+ s{(b+iA)(2^-i/4)‘^^+ z[i(z-i/2)A-bz](z^-i/4)^''^} + (?(|slb^) 

Equation (A23) is the classical first-order approximation for the flow past an 
ellipse of thickness b situated at angle of attack A relative to a uniform stream. 
It can be shown that the outer first-order solution satisfies the first-order flow 
tangency condition, Eq. (A15) , when the steady velocity is obtained from the outer 
approximation to the zeroth order solution, Eq. (A23) . It follows from Eqs. (A18) , 
(A23) and (A24) that the outer approximation to the pressure difference across the 
airfoil is 


Ap(x) 


-4a[(i/2-x) + b(i-3x)j(i/4-x^) 0 (|s|b^) (A25) 


Therefore, to lowest order Aplx)!^ for the displaced airfoil depends on airfoil 
thickness, b, but not on mean angle of attack, A . 

Near the leading edge of the airfoil (zj^g 0 the outer solution assumes 

the form 
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I - (A+ib/2) + ... . 


(A26) 


(i-b) 


dw 

dz 


0 


and 




-s(A + ib/2)z 


-3/2 


LE 


n- 0.1 


- 1/2 


LE 


(A27) 


2 

where = z + 1/2 = z^, + b /4. Thus, to lowest order the inner and outer solutions 
have the same functional form near the leading edge of the airfoil. The inner solu- 
tions depend on the variable Zq which is measured from the center of the leading 
edge circle, while the outer solutions depend on z^^g which is measured from the 
leading edge of the airfoil. 


Closure 

The foregoing model-problem study reveals several important features of pertur- 
bation solutions to airfoil type flows. In particular, the solution for the flow 
past a displaced airfoil can be represented as an asymptotic series in which the 
zeroth order solution describes the flow past the airfoil at its original location 
and the remaining terms account for the effects of the displacement. This series 
will converge to the exact solution for the flow past the displaced airfoil if the 
displacement is small compared to the edge radii of the airfoil. If this criterion 
is not met, but the displacement is still sufficiently small, the asymptotic series 
will converge everywhere except near the edges of the airfoil. In this case the 
perturbation approximation is singular. In either event the first-order term of 
the asymptotic series is a solution of a linear boundary value problem in which 
the airfoil surface boundary condition (i.e., the flow tangency condition) is satis- 
fied on the original airfoil surface. The zeroth-and first-order solutions are 
analytic for blunt-edged airfoils, but these potentials, or their derivatives to 
some order, are singular at sharp leading and trailing edges. Finally, near blade 
edges, the zeroth-and first-order solutions are essentially functions of a position 
vector emanating from the center of the edge circle. 
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Numerical Implications 


One purpose in conducting the model-problem study is to provide guidance for 
the numerical resolution of unsteady cascade flows. First-order results depicted 
in Figs. 21 through 25 indicate that with the exception of small neighborhoods 
surrounding the blade edges, the first-order unsteady flow could be captured on 
a relatively coarse rectilinear type mesh which covers the extended blade passage 
region; e.g., the mesh shown in Fig. 3. Although this mesh is generally well suited 
for approximating conditions on the boundaries of the extended blade passage region, 
and for capturing the unsteady flow throughout Tnost of the solution domain. It is not 
suitable for resolving the flow in the vicinity of blunt blade edges. As a result 
solutions to the unsteady-cascade boundary value problem were initially reported 
for sharp-edged profiles (Refs. 7 and 11). For a blunt-edged airfoil, a solution 
on the cascade mesh must be supplemented by local solutions on very fine meshes 
(e.g., Fig. 4) in regions where velocity gradients are large. The coarse mesh solu- 
tion is then used to provide conditions on the outer boundary of the local mesh. A 
prerequisite for this approach is the ability to determine a reasonably accurate 
unsteady solution on the coarse mesh. 

The present authors' early attempts to resolve the unsteady boundary value 
problem for blunt leading-edged, compressor-type, airfoils on the cascade mesh 
(Fig. 3) resulted in spurious unsteady pressure predictions over the entire airfoil 
including violent pressure oscillations near the leading edge. Results of the model 
problem study reveal the source of difficulty to be the partial inclusion of local 
phenomena in the coarse mesh calculation. Since this mesh should only be used to 
calculate an "outer" approximation to the unsteady flow field, local steady flow 
behavior should be excluded from the governing differential equation, Eqs. (9, 10), 
and surface boundary condition, Eq. (11), for the coarse mesh calculation. This is 
accomplished by not placing mesh lines too close to a blunt leading edge, and by 
replacing the local steady flow solution by an approximate outer solution in the 
vicinity of blade leading edges. That is, the functional dependence of steady 
potential derivatives on a position vector emanating from the center of the leading 
edge circle, is replaced by a functional dependence on a position vector emanat- 
ing from the leading edge, to produce singular, but monotonic (near = 0) 

behavior in steady potential derivatives along radial lines emanating from the lead- 
ing edge. Such behavior poses no difficulty to the calculation of coarse mesh solu- 
tions. The unsteady numerical solution on the coarse mesh is then akin to an outer 
unsteady solution, and as such, it provides reasonably accurate boundary condition 
information for determining an inner solution on a dense local mesh surrounding a 
blunt leading edge. Similar considerations should be applied to calculate the 
unsteady flow past a cascade of airfoils with blunt trailing edges. 
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Fig. 1 Two-Dimensional Oscillating Cascade with 
Finite Mean-Flow Deflection 
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Fig. 3 Cascade Mesh 
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Fig. 5 Unsteady Pressure Difference near the 
Blade Due to Bending 
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Fig. 6 Unsteady Pressure Difference near the Leading Edge of a NACA 0012 
Blade Due to Torsion 
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Fig. 8 Effect of Inlet Mach Number on Unsteady Response Due to Bending for an 
Unstaggered Cascade of NACA 0012 Airfoils 


47 




UNSTEADY PRESSU, .E DIFFERENCE COEFFICIENT, Ap 



Fig. 9 Effect of Inlet Mach Number on Unsteady Response Due to Torsion for an 
Unstaggered Cascade of NACA 0012 Airfoils 
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Fig. 11 Effect of Inlet Mach Number on Unsteady Response Due to Torsion for an 
Unstaggered Flat Plate Cascade 


50 





REAL UNSTEADY PRESSURE COEFFICIENT, Re {p| 


30 





Re |C^^| 

0° 

-6.46 

1.80 

5° 

-6.42 

1.81 

10° 

-6.35 

1 

1.85 

I 


0.2 

0.3 


SQUARE-ROOT DISTANCE FROM LEADING EDGE.vCT 


Fig. 12 Effect of Incidence on Real Pressure Distribution near the Leading 
Edge of a NACA 0012 Blade Due to Torsion 
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Fig. 13 Effect of Incidence on Imaginary Pressure Distribution near the 
Leading Edge of NACA 0012 Blade Due to Torsion 



































CO 


_] 

< 


O 

CL 

> 

t 

o 

o 

_] 

LU 

> 


LU 

o 

< 

ll 

q: 

3 

CO 



REAL COORDINATE, x 

UPPER SURFACE LOWER SURFACE H 

TE LE TE 


Fig. 21 First-Order Potential on Airfoil Surface for Incompressible Flow 
Past Flat Plate and Elliptic Airfoils 
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Fig. 22 First-Order Normai-Veiocity on Airfoil Surface for Incompressibie Flow Past 
Flat Plate and Elliptic Airfoils 
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Fig. 23 First-Order Normal-Velocity on Airfoil Surface near Leading Edge 
for Incompressible Flow Past Flat Plate and Elliptic Airfoils 




PRESSliRE DIFFERENCE, Ap 



Fig. 24 First-Order Pressure Difference for incompressibie 
Flow Past Flat Plate and Elliptic Airfoils 








Fig. 25 First-Order Pressure Difference near the Leading Edge for 
Incompressible Flow Past Flat Plate and Elliptic Airfoils 
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